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Abstract 

It is shown that the interacting Pais-Uhlenbeck oscillator necessarily leads to a descrip- 
tion with a Hamiltonian that contains positive and negative energies associated with two 
oscillators. Descriptions with a positive definite Hamiltonians, considered by some authors, 
can hold only for a free Pais-Uhlenbeck oscillator. We demonstrate that the solutions of 
a self-interacting Pais-Uhlenbeck oscillator are stable on islands in the parameter space, 
as already observed in the literature. If we slightly modify the system, by considering a 
sine interaction term, and/or by taking unequal masses of the two oscillators, then the sys- 
tem is stable on the continents that extends from zero to infinity in the parameter space. 
Therefore, the Pais-Uhlenbeck oscillator is quite acceptable physical system. 

1 Introduction 

The ordinary gravity, described by the Einstein-Hilbert action (containing the cur- 
vature scalar R), is not renormalizable. The higher derivative gravity, with R + R 2 , 
is renormalizable. But higher derivative theories are considered as problematic, be- 
cause according the the Ostrogradski formalism [lj they contain negative energies 
which, according to the wide spread belief, automatically imply instabilities at the 
classical and quantum level. It is often stated that at the quantum level such a theory 
implies negative probabilities due to ghosts, and is therefore not unitary. Whether 
a theory implies negative probabilities and positive energies, or vice versa, positive 
probabilities and negative energies, depends on choice of vacuum and corresponding 
creation and annihilations operators [2] -[5]. 

A model for a higher derivative theory is the Pais-Uhlenbeck oscillator [7J . It 
has been studied by many authors, because understanding the issues concerning its 
stability, could pave the way towards quantum gravity. Smilga [S] has found that there 
are islands of stability of the classical solutions of the interacting Pais-Uhlenbck (PU) 
oscillator. Recently, Mostafazadeh [S] has found a Hamiltonian formulation of the PU 
oscillator that yields a stable and unitary quantum system. Other authors [rO]-[To] 
have also arrived at the positive definite Hamiltonians for the PU oscillator. In this 
paper we show that such a description holds for a free PU oscillator only, but not for 
a self-interacting one. In the latter case one has to describe the PU oscillator by the 
second order Lagrangian and employ the Ostrogradski formalism. Then, as it is well 
known, the PU oscillator can be written as a system of two coupled oscillators, one 
with positive and the other one with negative energy. Stability issues arise, because 



the energy can flow from one to the other oscillator in such a way that their kinetic 
energies escape into positive and negative infinity, respectively, while the total energy 
of the system remains constant and finite. 

In this paper we first show explicitly that, in general, a self-interactin PU oscillator 
cannot be described by a positive definite Hamiltonian. Then we consider numerical 
solutions of the PU oscillator with the quartic self-interaction. Unless the coupling 
constant or the initial velocity is too high, the (classical) system is stable. So we 
indeed have islands of stability as observed by Smilga [8], and recently by [17]. Then 
we consider two modifications of the PU oscillator, (i) We replace the quartic inter- 
action term with a term that contains the forth power of sine. We show numerically 
and analytically that such a modification gives continents of stability, (ii) Instead of 
taking equal masses of the two oscillators, we consider the case in which the masses 
are different. Such a modified system is in fact just like the PU oscillators, only the 
coefficients in front of the terms are changed, and a non-linear term is added. Again 
we obtain stability for the vast range of parameters. Moreover, regardless of how 
high the initial velocity is, the solution is stable. Such behavior of classical solution 
implies that the quantum system is stable as well [17J. 

2 The Pais-Uhlenbeck oscillator as a system of two 
oscillators 

The Pais-Uhlenbeck oscillator is given by the 4th order differential equation 



As observed by Mostafazadeh (9] , this can be written as the system of two oscillator: 





which gives 




x + — piy = 



(3) 



y + ^y- p2% = o 



(4) 



where pi, p 2 , Pi, P2 are real constants. 




Inserting this into Eq. (|4]), we obtain 




(5) 



Comparison of the latter equation with ([2]) gives the relations 

Pi + P2 = w\ + U 2 



(6) 



1 We use here a different notation for coefficients. 
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PlP2 ~ Plp2 = (7) 

The solutions is 



.2 

1,2 



§(//! + /i 2 ) ± W(^+ ~ 4(/Xl/X 2 - PlP2)- 



Let us now find possible Lagrangians corresponding to the equations of motion 

(§,©■ 
Case I. 

Assuming the Lagrangian 

L = 1(A 2 + y 2 ) - ^ lX 2 + fi 2 y 2 - 2 Pl xy), (9) 
we obtain the equations of motion (|3]),(|4j), if 

P2 = Pi- (10) 

Then from Eq. (|8| we have 



U l = + ^ 2 ) ± I. - ^2 + 4 p 2 . (11) 



We see that for a large range of the coefficients /ii, /i 2 , pi, the squared frequencies uf 
and bj\ are positive. Then co>i, u 2 are real, in which case we have oscillating motion. 
The Hamiltonian is 

H= l(p 2 x + p 2 y ) + ±(fnx 2 + fi 2 y 2 - 2 Pl xy), (12) 

where p x = dL/dx = x, and p y = dL/dy = y. By performing a rotation in the 
(x, y)-space, 

x = x cos a + y sin a 

y' = — x sin a + y cos a (13) 
with the accompanying rotation of momenta, 

Px' = Px cos a + p y sin a 

Py' = ~Px srn a + Py cos a ; (14) 



the Hamiltonian (12) can be diagonalized. By comparing the new Hamiltonian 

H=\{p 2 x ,+p 2 y ,) + \{ax> 2 + by> 2 ), (15) 



with the old one (12), we obtain the system of three equations for the unknowns a, 
6, a: 



a cos 2 a + b sin 2 a = fii 

a sin 2 a + b cos 2 a = fx 2 (16) 
(a — 6) cos a sin a = p\. 



The solution is 

a = !(//! + n 2 ) + \\J Oi - P2) 2 + 4pf = u\. (17) 



b = §fai + (Mi) - \\J^i - /i 2 ) 2 + 4p? = u; 2 . (18) 

cos2a=-=^L ~ /X2 (19) 
V0"i " /^) 2 + 4p? 

The a, 6 are just equal to the squared frequencies w 2 , of the PU oscillator. This 
can be directly verified by inserting the expressions (16) into the equations of motion 
([5]) and using (10). 

In the new coordinates, the system is described by the Lagrangiarj^] 

L = \{x> 2 + y' 2 ) - \{ulx' 2 + coly'% (20) 

and the Hamiltonian 

H=\(x' 2 + y' 2 ) + \(uj 2 1 x' 2 + u J ly' 2 ). (21) 

The energy of this system is positive. It is remarkable that when we diagonalize the L 
and H for a system of two oscillators ^ and Q, we obtain two different frequencies, 
U2 and 0J2, that correspond to those occurring in the PU oscillators. 
Case II. 

Alternatively, we may assume that the Lagrangian is 

L = \{x 2 - y 2 ) - \{^x 2 - fx 2 y 2 - 2 Pl xy), (22) 
This gives the equations of motion ^ and ^ if 

P2 = -pi ■ (23) 

Inserting this into Eq. (pi), we obtain 



to 



I = 1^ + M2 ) zp 1 _ - Ap\. (24) 



The frequencies ui, U2 are real if (/xi — /X2) 2 > 4p 2 and fix + /X2 > a/ (/ii — /X2) 2 — 4pf . 
The Hamiltonian is 

# = I fa* - pj) + \(Pix 2 - A^ 2 - 2p!xy) . (25) 
By performing the hyperbolic rotation in the (x, y)-space, 

x' = x cosh a + y sinh a 

y = x sinh a + y cosh a (26) 



2 The author [5] has also arrived at such a system of two uncoupled oscillators straightforwardly 
from Eq. ^ , by using a different chain of substitutions of variables. See also the procedure of 
refs. ITS]. 
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with the accompanying rotation of momenta, 



Px' 
Py' 



p x cosh a + p y sin a 
p x sinh a + p y cosh a, 



(27) 



the Lagrangian (22) and the Hamiltonian (25) become 

u/2 



L - 
H 



[x ■ 



2 (x 



:'2\ 



V 



( 2 12 
(OJiX 



2„.i2\ 



u 2 y 



:'2\ 



(28) 

!/-) + HU'l-r- - ^o//'-J. (29) 

Again, the diagonalized Lagrangian and Hamiltonian contain the frequencies lu±, uj 2 
of the PU oscillator. 

Now we have the relations 



u\ cosh 2 a — ui^, sm h 2 a = fii 
—cuf sinh 2 q. + Cc>2 cosh 2 a 



w 2 ) cosh a sinh a 



M2 

-Pi- 



(30) 
(31) 
(32) 



The energy of the system is either positive or negative, depending on which degree 
of freedom is more excited. 

Cases I and II show that the PU oscillator can be described as a system of two 
oscillators whose Hamiltonian is either (21) or (29). Case I means positive definite 
Hamiltonian, whereas Case II means indefinite Hamiltonian. 



3 Self-interacting PU oscillator 

3.1 Equations of motion and the Lagrangian 

We have seen that the PU oscillator can be described as a system of two oscillators 
Q and Q that can be written in the explicit uncoupled form 

x' + cofx' = (33) 

y' + u 2 2 y' = (34) 

For real u)\, this is an oscillating system, regardless of whether for the corre- 
sponding Lagrangian we take (20) or (28). Both Lagrangians are equally good for 
describing the PU oscillator JTUJ H~5] . 

If we include an interaction between the x' and y', then energy can be transfered 
between those two degrees of freedom. Then it does matter which Lagrangian we 
take. 

(i) Let us first consider the following Lagrangian that is an extension of (20) 
(Case I): 



\{x' 2 + y' 2 ) - \{u\x' 2 + u 2 2 y' 2 ) - ±(x> + ,V (30) 



The corresponding equations of motions are 

x' + u\x' + \(x' + y') 3 = (36) 

y' + u 2 2 y' + X(x' + y'f = (37) 
Introducing the new coordinates 

x' +y' x' — y' 

u = — 7 =- , v = — (38) 

y/2 ^2 



we have 



where 



L = \{u 2 + v 2 ) - + u 2 2 )(u 2 + v 2 ) + 2K 2 - u 2 2 )uv] - ^u 4 (39) 

u + fxiu - piv + Am 3 = (40) 

v + p 2 v - p\u = 0, (41) 

A*i = ^2 = + ul) , -pi = -Wa)- ( 42 ) 



Eliminating it, we obtain the 4th order differential equation for v. 

« (4) + (Mi + + - p\)v + Xpi(v + p 2 vf = 0, (43) 

which is just that of the PU oscillator with an extra non linear term. 
Similarly, by eliminating v, we obtain 

d 2 

M (4) + (/X! + p 2 )ll + (/X!/X 2 - p 2 )iz + /i 2 Ali 3 + A— (u 3 ) = 0, (44) 
which is also the PU oscillator with a non-linear term. 



(ii) Let us now consider the Lagrangian that is an extension of (28) (Case II): 



L = \(x> 2 - y' 2 ) - \{uj 2 iX ' 2 - oo 2 y> 2 ) - ~(x' + y'f (45) 
The corresponding equations of motions are now 

x' + ulx' + \(x' + y'f = (46) 

y> + UJ ly>_ X {x' + y'Y = U (47) 

Notice the minus sign in the second equation. 

In the new variables it, v, defined in Eq. (38), we have 

L = iiv- - ujI){u 2 + v 2 ) + 2(tu 2 + uj 2 2 )uv] - ^u 4 (48) 



6 



u + HiU — p\V = 
v + \i 2 v — p\U + Am 3 = 0, 



where \i 2 and pi are given in Eqs. (42). By eliminating v, we obtain 
m (4) + (/ii + fi 2 )u + (p>ifi 2 - p\)u + pi Am 3 = 0. 



Using (42), the latter equation reads 



it: 



+ {u\ + u 2 )u + uj\uj\u — 



(w? - col)\u 3 = 0. 



(49) 
(50) 

(51) 
(52) 



Now we obtain the equation of motion for the PU oscillator with a self-interaction 
term. The second order Lagrangian that gives the fourth order equation of motion 
( 52 ) is that of the PU oscillator with a quartic self-interaction term: 

L = \ [u 2 - (col + cuj)u 2 + uxulu 2 ] +\{ul- uj 2 2 )\u a . (53) 



Notice that if u\ = u?, then p\ = 0, and v in (49) cannot be expressed in terms of u. 
Consequently, in the case U\ = 002, the system (49), (50) does not give the equation 
(51) for the PU oscillator. 

The Ostrogradski second order formalism then leads us to the phase space La- 
grangian [lU [IH1 EE] that is equivalent to (53), 

L = p u ii + p q q - H, 

where 



(54) 



H = p u q + l [p 2 q + (w? + tu 2 2 )q 2 - uju 2 2 u 2 ] + \Au\ (55) 
Here we introduced A = \{ui 2 — w|)A. The latter Hamiltonian can be transformed 



into 



H' = \{pl - pi) + H^x' 2 - ujy' 2 ) + X(x' + y')\ 



The phase space Lagrangian 



L = p x ix +p y >y - H 



(56) 



(57) 



is equivalent to (45). This can be directly seen by using the equations of motion 
p x i = x' and p y i = y', and eliminating p x >, p y > from (57). 

Therefore, the correct procedure is to start from the 4th order self-interacting 
Lagrangian (53) and to employ the Ostrogradski formalism. The Hamiltonian so 
obtained has positive and negative values. The procedures discussed in Refs. [HI US], 
and also in Refs. [10] [IS] , have limited validity, because they do not consider an 
interaction term. They are valid descriptions of the PU oscillator in the absence of 
an interaction, but not if one switches on an interaction. 
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3.2 Solutions 



The Lagrangian of the form (45) is usually considered as unsuitable for for physics, 
because it implies indefinite Hamiltonian, with positive and negative energy states. 
The interaction term that mixes the two types of states leads to instabilities. But as 
pointed out in Refs. [HIE], there exist islands of stabilities. We show this explicitly by 
solving numerically the equations of motions (46), (47). In Fig. 1 there are examples 
of such calculations, done by MATHEMATICA. In all examples we take uf — 1 and 



1.5 



y 











Figure 1: Solutions of Eqs. (46 1 (47 1 for different values of the coupling constant A and different 
initial conditions. Left and middle: the trajectories in the (x',y') space. Right: The kinetic energy 
x' 2 /2 as function of time. The oscillations within the envelope are so fine that they fill the diagram. 



We see that the system is stable for sufficiently small coupling constant A and the 
initial velocity x'(0), y'(0). If A is too high, the system is unstable (Fig. 2, up). 

Similarly, the system is unstable at too high velocities (Fig. 2, down). Close to the 
critical value of A, the system seems to be stable for long time, but then it escapes 
into infinity (Fig. 3). A similar behaviour occurs close to the critical value of the 
initial velocity. 

By just slightly decreasing the coupling constant from that of figure 3, A = 
0.02299, to A = 0.0229, the system appears to be stable. We checked its stabil- 
ity up to t = 2664, but we do not plot the solutions here, in order to not crowd the 
paper with too many figures. 
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Figure 2: Up: By increasing the A, the system becomes unstable. The trajectory and the kinetic 
energy escape into infinity. Down: Similarly, by increasing the initial velocity, the system also 
becomes unstable. 




Figure 3: At certain values of A and the initial conditions, the system behaves stably for a long time, 
before it finally escapes to infinity. The total energy E to t remains constant within the numerical 
error. 



The interaction potential j(x'+y') runs into infinity. More realistically, it should 
not run into infinity, but there should be a cutoff. As a more realistic coupling term 
let us consider |sin 4 (x' + y'), that leads to the Lagrangian (53 ) in which w 4 is replaced 
by sin 4 w. The equations of motion are then 



H + u{x' + Asm 5 (x' + ?/')] cos (x' + y') = 0, 



y + u 2 y — Asin (x + y ) cos (x + y 







(58) 
(59) 



Such system is stable at all values of A > and initial velocity. We have checked this 
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by performing many numerical runs. In (Fig. 4) we give two examples of numerical 
solutions. Later we will demonstrate also analytically why the solutions of the system 



(58), (59) are stable 



A=0.22 

x'(0)=0, y'(0)=l 
x'(0)=l, jf(0)=0 






A=l 

x'(0)=0.3, y'(0)=l 
/2 x'(0)=l, y'(0)=-0.5 



- t 



Figure 4: Solutions to the equations of motion ( 58 1 , ( 59 1 , in which the quartic interaction h (x'+y') 4 



is replaced by jsin (x' + y'). The system is now stable for all positive values of A. 



Another possible generalization is in replacing the Lagrangian (|45|) with 

•/2\ 



L = \{rn\x a — m 2 y' 



A 



iKV 2 -c 2 V 2 )-^(x' + 2/ ') 4 



(60) 



where mi and m>2 are now two different "masses". In terms of the variables u, v, we 
have 



L = \ [m(u 2 + v 2 ) + 2Miiv + p x (u 2 + v 2 ) - 2/i X w] - 



A 



-u 



where 



m = \{jn\ — m 2 ), M = \{fn\ + m 2 ) 
The equations of motion are now 

mil + Mi) — p\u + \i\V + Am 3 = 

mv + Mil — p\v + \i\V = 0. 
The corresponding 4th order equation is 



(61) 
(62) 

(63) 
(64) 



u 



d 2 

(4) M(M 2 - m 2 ) + u2M(fnM + p^) + uM(p\ - p\) + Mp 1 Xu 3 - MmX— (m 3 ) = 

(65) 
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This is a deformed version of the equation (51) for the interacting PU oscillator. By 
taking m = 0, M = 1, we obtain the ordinary PU oscillator of Eq. (51 ). 
Examples of numerical solutions to the equations of motion 



m\x' + uj x x' + \{x + y) 



A3 



m 2 y + u 2 y' - \{x' + y') 



/\3 



0. 



(66) 
(67) 



derived from the Lagrangian (60), are given in Fig. 5. Whilst in the case of equal 
masses, mi = m2 = 1, the system is unstable at A = 0.03 and higher, we see that 
for different masses, mi < 1712, the system is stable regardless of the values of A > 
and initial velocities. This has been confirmed in many numerical runs that we have 
done. Only a small sample is shown in Fig. 5. That for unequal masses the system 
becomes stable we previously observed in Ref. [20], where we studied an analogous 



system of two oscillators, but with a different coupling term, namely, 



x' 



2\2 



-y 2 ) 



which is a special case of that considered in Ref. [TTj . However, such a coupling term 
does not correspond to a quartic self-interaction of the PU oscillator. 



A=5, mi=0.7, T7i2=1.3 
s'(0)=0.3, y'(0)=l 
x'(0)=4, j/'(0)=-0.5 
wi=l, o>2=VX5 




A=500, mi =0.7, m 2 = 1.3 
x'(0)=0.3, y'(0)=l 
i'(0)=4, y'(0)=-0.5 
,„ wi=l, oj 2 =\/T5 



A=5, rrti=0.7, m2=l-3 
x'(0)=0.3, 2/'(0)=l 
i'(0)=40, y'(0)=55 

U3l=l, OJ 2 =\/l.5 



A=5, mi=0.99, m 2 =1.01 
x'(0)=0.3, y'(0)=l 
x'(0)=l, y'(0)=-0.5 
a)i=l, aj2=\/l^5 



A=5, mi=0.7, m 2 = 1.3 
x'(0)=0.3, y'(0)=l 
i'(0)=40, y'(0)=55 

U)l =0)2 =0 

Ifl|l||Nl|l|l|ln]|l|l|M 



A=5, mi=0.7, m2=1.3 
x'(0)=0.3, y'(0)=l 
x'(0)=4, y'(0)=-0.5 

O)l=O) 2 =0 



II II 



"t 



I ..'.'"'j .j. 1 '' 1 .' .i 1 " 1 '. J'.'.', i. 1 '! 1 ' , 



t 



t 



Figure 5: Solutions of Eqs. (66 1 (67 1 for different values of the coupling constant A and different 
initial conditions. We show here the kinetic energy i' 2 /2 as function of time. 

To see why the system with different masses is stable, let us inspect the equations 



of motion (66), (67). We already know that at small values of A, the system is stable. 
At large values of A, we can neglect the terms u\x' and io\y' '. Equations of motion 
are then 

x'+ — \{x' + y') 3 = 0, (68) 
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y' - —X(x' + y'f = 0. (69) 
Taking the sum and the difference of the latter equations, we obtain 

i+(- -W = 0, (70) 



1 


1 


nil 


m 2 


1 


1 

f- — 


mi 


m 2 



V+[— + —)\e = 0, (71) 



where £ = x' + y' and r\ = x' — y' . 



In the case of unequal masses, nil < m 2 , A > 0, Eq. (70) describes the quartic 



oscillator with the quartic potential, |£ 4 , which has stable, oscillatory solutions. 
Then, Eq, ( 71 ) also has stable, oscillatory, solutions. Stability is maintained in the 
presence of the terms ufx' and u 2 y'. 

In the case of equal masses, mi = m 2 , Eq. (70) becomes £ = 0, with the solution 
£ = £o + c\t. Then the general solution of (71) is a runaway function 

V = - — ^o + cit) 5 + c 2 t. (72) 
mi 20 

In the presence of the terms ufx and u 2 y, the above runaway behavior is modulated 
by oscillations. 

Solutions to the system described by the Lagrangian (60) are stable, if mi > m 2 , 
A > 0. If mi = rn2, then the solutions are stable at sufficiently small A, whereas at 
higher values of A, they are unstable. (Fig. 2). 

If instead of \(x' + y') we take the interaction term ^sin 4 (x' + y'), we have stability 
even in the case m x = m 2 . The equations of motion are then 

2 

£ = 0, n+ — Asin 3 £cos£ = 0, (73) 
mi 

the general solution being 

£ = £o + cit, ») = -— -j— Asin 4 (£ + cit) , (74) 

mi 4ci 



V = ~ — ^2 X [ 12 (& + cit) - 8 sin (2(6, + cit) + sin (4(£ + erf)] (75) 



2 1 

m~il28cj' 

This solution is stable in the sense that the velocity and the kinetic energy remain 
finite. The coordinates £, rj, or equivalently, x', y' proceed with time, on average 
linearly, into infinity. The velocity thus oscillates around a constant velocitjj^] If we 
include into the potential also the terms \(J\x' 2 and ^ufy' 2 , then the coordinates do 
not escape into infinity, but they oscillate. 



3 In the quantized theory, to such modulated uniform motion there corresponds a modulated 
traveling wave, or uniformly moving wave packet. 
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4 Discussion 



It has been shown by some authors [9j[T6] (see also [10] [IS]) that the Pais-Uhlenbeck 
oscillator can be described as a system of two degrees of freedom with positive definite 
Hamiltonian. We point out that this holds for the free PU oscillator only and that 
one cannot include a couppling term such that the system would be equivalent to the 
PU oscillator with a quartic or similar self-interaction term. The interacting Pais- 
Uhlenbeck oscillator must be described, as usually, by the second order Lagrangian. 
The Ostrogradski formalism then leads to the indefinite Hamiltonian, with positive 
and negative energies. An equivalent system is that of two oscillators described by 
the equations of motion (46), (47), derived from the Lagrangian (45). We have studied 
numerical solutions to the latter system for various coupling constants A and initial 
velocities. Solutions are stable below a critical value of A and initial velocity. We then 
considered two modifications of the Lagrangian that drastically increase the range of 
stability. 

Firstly, we replaced the quartic interaction term | (x'+y') 4 , that runs into infinity, 
with the term j sin 4 (x' + y') that is finite for all x', y' . Then, instead of islands of 
stability, we obtain the continent of stability. Fig. 4 shows that now the system is 
stable even at A = 5, whereas with the quartic interaction it was unstable already 
at A = 0.03. We have done many numerical runs with higher values of A, even 
with A = 500, and the solutions were always stable. By inspecting the equations of 
motion, we also found analytically that such interacting system is indeed stable for 
any positive A, and for any initial velocity. 

Secondly, we replaced the kinetic term \{x' 2 — y' 2 ) with \{ra\x' 2 — m 2 ?/ 2 ), and 
considered the case in which the "masses" mi and m 2 are different. If m\ < rri2, 
A > 0, and uf < u> 2 , the system is stable for all finite positive values of A and for 
all finite positive or negative initial velocities x'(0), y'(0). Analogously, the system is 
stable if mi > m 2 , A < 0, and u\ > ui 2 . 

Our findings invalidate the generally held belief that the Pais-Uhlenbeck oscillator 
in the presence of an interaction is unstable, and therefore problematic. There are 
vast regimes of stability that hold for all initial velocities. This has consequences 
for the quantum PU oscillator. Namely, stability of a classical system does not 
necessarily imply stability of the corresponding quantum system, because the latter 
system can tunnel through a potential barrier and then role down the potential. But 
if a classical system remains stable, regardless of how high is the initial velocity, then 
also the quantum system is stable. We conclude that the Pais-Uhlenbeck oscillator 
with a suitable self-interaction is quite acceptable from the physical point of view. 
Since the PU oscillator is a toy model for higher derivative gravity, we expect that 
also the negative energy problems of the latter theory could be resolved along similar 
lines as investigated in this paper. 
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